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Abstract 

In this lecture I present some of the new developments concerning the use of 
Pade Approximants (PA's) for resumming perturbative series in QCD. It is shown 
that PA's tend to reduce the renormalization scale and scheme dependence as 
compared to truncated series. In particular it is proven that in the limit where the 
(3 function is dominated by the 1-loop contribution, there is an exact symmetry 
that guarantees invariance of diagonal PA's under changing the renormalization 
scale. In addition it is shown that in the large (3q approximation diagonal PA's can 
be interpreted as a systematic method for approximating the flow of momentum 
in Feynman diagrams. This corresponds to a new multiple scale generalization of 
the Brodsky-Lepage-Mackenzie (BLM) method to higher orders. I illustrate the 
method with the Bjorken sum rule and the vacuum polarization function. 



1 A lecture given at the Cracow School of Theoretical Physics, Zakopane, May 30 - 
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I will talk about resummation of perturbative series in QCD The basic 
question I deal with is how to use finite order perturbative calculations in 
QCD to make unambiguous theoretical predictions, with controlled errors. 
As experiments improve one requires the theoretical predictions to be more 
accurate. However, in QCD it is very hard to get accurate predictions, ba- 
sically because the coupling constant is large. This leads to non-negligible 
non-perturbative effects as well as a badly divergent and renormalization 
scheme dependent perturbative series. 

In this talk I will show that PA's which start out as an alternative to 
a finite order perturbative series having the same formal accuracy, actually 
have an important advantage over the finite order series^. Through the re- 
summation of certain all-order effects related to the running of the coupling- 
constant, PA's become independent of the choice of the renormalization scale 
and therefore lead to more accurate and more reliable predictions. The ma- 
terial presented in this lecture appears in greater detail in ref. J|, |J. 

The outline of the talk is as follows: I will start by introducing the PA's 
method and the problem of renormalization scale dependence in QCD. I will 
mention some of the other ideas that were raised to confront the problem 
of renormalization scale dependence and show how PA's solve it in a most 
elegant way. Then I will address the question of what higher order effects are 
summed-up by PA. I will show that there is a direct interpretation of PA's 
in terms of approximating the momentum distribution of virtual gluons in 
Feynman diagrams. 

First, what are PA and how do I use them? I start with an effective 
charge related to some physical observable, written as a power series in x: 

S n = x (l + T\x + r 2 x 2 + • • • + r n x n ) (1) 

where x = a s /7r. A PA is constructed by writing a ratio of two polyno- 
mials such that when expanded back to a Taylor series, it gives the known 
coefficients r\ through r n : 

1 + CL\X + ... + (InX N 

1 + bix + ... + b M x ! 



P x[ n/m] =x - , ^ , 7T - ZM , N + M = n 



2 For a recent review on the nature of perturbative series and resummation techniques 
see ref. 

3 PA have various successful application in physics. Examples of applications to statis- 
tical physics and quantum field theory are listed in ref. [Q. Applications to QCD appear 
in ref. fll 
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There is a theorem that for any degree N in the numerator and M in the 
denominator such that N + M = n there is a unique PA function || . I will 
mainly deal here with diagonal PA's which are written, in my notation, as 
x[N/N + 1], having one power of x out of the brackets. I will soon come back 
to discuss PA's. 

A renormalized perturbative series in QCD is not expected to give exact 
predictions for measurable quantities due a few limitations. First, the series 
is divergent and not even Borel-summable. The resulting ambiguity is related 
to the existence of non-perturbative effects. Second, at any given order, the 
partial-sum depends on non-physical parameters, such as the renormalization 
scale. This also makes the prediction ambiguous. 

Let us concentrate on the renormalization scale dependence. Usually, 
when we calculate some observable R that depends on one external mo- 
mentum Q 2 in perturbation theory, we choose as an expansion parameter 
the renormalized coupling-constant at the external scale Q 2 . This "natu- 
ral" choice of fi 2 = Q 2 is, however, quite arbitrary. We can, just as well, 
use some other expansion parameter y = a s (e t Q 2 )/7c, where t ^ 0. The 
renormalization group equation 

-J = f3 x 2 + fax 3 + f3 2 x 4 + ■■■ 
at 

determines how the two couplants are related, 

x = y + (3 t y 2 + (fit 2 + p x t) y 3 + (> 3 t 3 + ^ ftt 2 + ftt) y 4 + ■ ■ ■ (2) 

and thus how the finite order series can be written in terms of y: 

S n {t)=y (l + r x y + f 2 y 2 + f 3 y 3 ...+ r n y n ) . (3) 

The new coefficients fj are different from the original coefficients r^, so as 
to compensate for the scale shift, such that the total effect is some residual 
dependence on t which is of the next, uncalculated order. Still, in QCD, since 
the coupling constant is large, the numerical difference due to the change of 
scale can be quite large. This limits the predictive power of the theory. 
Beyond two-loops, there is also the question of scheme dependence which 
can be parameterized by the higher-order coefficients of the (3 function, 
03 and onward. 
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In order to test the significance of this scale and scheme dependence, we 
studied the polarized Bjorken Sum-Rule. In fig. 1 the Bjorken effective 
charge at NNLO for Q 2 = 20 GeV 2 is plotted as a function of the renormal- 
ization group non-physical parameters: the coupling x = a s (fi 2 )/ii and the 
second coefficient of the /3 function: C 2 = fa/Po- We see that the surface is 
far from being flat. 

The same surface is drawn again in fig. 2, but here - as a contour plot. 
The thick lines are contours of equal effective charge. Large renormaliza- 
tion scheme dependence corresponds to large higher-order corrections, since 
these are required to compensate for the scale dependence. This observa- 
tion makes it clear that we should carefully choose the renormalization scale 
and scheme that we are using. In fig. 2 one can identify a region of rela- 
tively low renormalization scale and scheme dependence. Specific scales and 
schemes are chosen according to different criterions such as the method of 
Effective Charges 0, the Principal of Minimal Sensitivity M and the BLM 
scale-setting method |§. For the Bjorken Sum- Rule example (fig. 2), all of 
the above are located in the central region of low renormalization scale and 
scheme dependence. Note that in this case, MS , with /i 2 = Q 2 is not a good 
choice. 

Let's go back to eq. (0) that describes the scale transformation relating 
the coupling- const ants x and y defined at two different scales. If we assume 
that the 1-loop coefficient of the (3 function, j3 Q , is large enough, i.e. 

Po > PiX { 

for any i > 1, we can approximate the full relation by one that includes only 
the leading terms in (3q\ 

x ~y + (3 ty 2 + Pit 2 y 3 + /? V + • • • 

This can be written in a closed form: 

= y 
l-Poty 

It is important to realize that in the physical case of QCD with 3 to 5 
flavors, this approximation is good. Fig. 3 shows the renormalization scale 
transformation itself, namely the running coupling constant as a function of 
the scale. The dashed line is the best we know of the running coupling in 
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QCD (it includes the 4- loop effects), and the solid line is the 1-loop running 
coupling. The two are quite close and I shall use here the 1-loop formula. 

I now get to the main point, which is the independence of PA on the 
renormalization scale. We saw that partial-sums as usually written in per- 
turbation theory always yield different results in different renormalization 
scales: 

S n (0) ± S n (t) 

where 5„(0) refers to fi 2 = Q 2 as in eq. ([!]), S n (t) refers to fi 2 = e t Q 2 as in 
eq. (0), and x — yj (1 — Poty)- However, if we construct a diagonal PA from 
the series in x (eq. ([[])), 

, s 1 + d\X + ... + a N x N 

P.WN + 1]{X) - X 1+biX+ ^ +bN+iX N + l 

and independently, another PA from the series in y (eq. (0)), 

5 ( v 1 + aiy + ... + a N y N 

Py[N/N+i]{y) = y ~ — — ~ — 

1 + hy + ... + b N+1 y N+L 

we will get the same result in both: 

Px{N/N+l](%) = Py[N/N+l](y)- 

This is due to the mathematical property of diagonal PA's: they are invariant 
under homographic transformations of the PA argument (x — > x/(l + Kx), 
see [|, Q). We know that the all-order result does not depend on the renor- 
malization scale. The fact that diagonal PA are invariant suggests that they 
correctly resum certain all-order effects that are related to the running of the 
coupling-constant. 

Non-diagonal PA are not exactly invariant. However, on the global level 
(for large scale shifts t) they always have a reduced scale dependence [|]. 
Going back to the example we examined above, namely the NNLO Bjorken 
sum-rule, we show in fig. 4 the a; [0/2] PA. Clearly (compare with the partial- 
sum of fig. 1) the renormalization scale and scheme dependence is almost 
completely eliminated! 

Non-diagonal PA's may be dangerous, since specific renormalization scales 
and schemes are sometimes particularly deviant, as in the example of the 
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x[l/l] PA for the Bjorken sum-rule shown in fig. 5. Therefore it is best to 
use a diagonal x[N — 1/N] PA. 

I now consider the question of what higher-order contributions are summed- 
up by diagonal PA's. It turns out that we can get some rigorous results [|J 



if we limit ourselves to the "large /?o" approximation (T^, IT[ |Tj|, where 
only the leading term in /3q in each perturbative coefficient is taken into 
account. This approximation corresponds to summing certain higher-order 
contributions that are related to the exchange of one virtual gluon. I use 
Neubert's formulation [10|, where resummation is achieved by using the run- 



ning coupling-constant at the vertices. The resummation integral is then a 
weighted average of the coupling-constant at all scales: 

w(k)a s (k 2 )d 4 k = / p(s)x v (e s Q 2 )ds 

J —oo 

where w(k) is the Feynman integrand and s = In {j^s)- The superscript V 
stands for the V-scheme which is the most convenient renormalization scheme 
for my purposes. While a specific scheme is used here in order to simplify 
the formulae, it is important to understand that the above resummation 



integral is scheme-invariant [10|. The function p(s) describes the distribution 
of momentum of the exchanged gloun, and x v (e s Q 2 ) describes the interaction 
strength as a function of the momentum. Using a 1-loop formula for x v (e s Q 2 ) 
I get: 

Clearly the integral includes contributions from an infinite set of diagrams. 
The exact distribution function (in the large /3q approximation) has been 
calculated for a few observables, such as the vacuum-polarization D-function 



[ |13| , I0| which I shall use here as an example. A representative diagram 



is the following: 
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In this particular diagram the exchanged gluon is dressed by fermion 
loops. However, gluonic corrections that are related to the 1-loop running of 
the coupling are resummed in the above technique just as well. 

Of course, the resummation integral is not well defined, due to the inte- 
gration over the infra red pole in the 1-loop formula for the running coupling 
(Landau pole). This is how infra red renormalons appear in this formulation. 
I will not deal here with the renormalon ambiguity which cannot be settled 
completely within perturbative QCD, but rather use the resummation inte- 
gral to study the PA method. The general methodology is to assume that 
only a first few coefficients in the perturbative series are known, construct 
a PA basing on this limited information and then compare the PA with the 
exact all-order result. We will see below that this comparison can be done 
also on the level of the momentum distribution function, since PA's can be 
interpreted as what one obtains by replacing the continuous momentum dis- 
tribution function with a particular discrete distribution. 

It was found empirically that the momentum distribution function in the 
large /3q approximation is a non-negative function in many physical examples 
[ PH , §]. This justifies a posteriori the probabilistic interpretation implied by 
the name 'momentum distribution'. If indeed p(s) > for any s, then the 
resummation integral defines a so-called Hamburger function ||: 

/OO y /*00 y 

p(s)-—ds = / — - d<f>(s). (4) 
-oo 1 + SZ J —oo 1 + SZ 

where z = /3oX V (Q 2 ) and </>(s) is the indefinite integral of p(s). The pertur- 
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bative coefficients are moments of the distribution function: 

/oo 
-oo 

for i > 0. 

There is a theorem || that guarantees that for a Hamburger function, a 
z[N — l/N] PA constructed from the partial-sum: 

2N-1 

* E fit-*)* 

8=0 

can be written as: 

f(z) ~z[N-l/N] =ETT^ 

i=i 1 "I" </i^ 

with g« rea/ and r; > for i = 1, 2, ...iV. 

Through this decomposition of the PA function, together with eq. (f|) one 
realizes that the PA corresponds to approximating the all-order continuous 
distribution function by a sum of iV weighted (^-functions: 

N 

Pn(s) = E r ^( s ~ 

i=l 

or, equivalently, its indefinite integral <fi(s), by a piece- wise constant function 
composed of N steps: 

N 

<I>n(s) = E r ^( s ~ Qi)- 

i=l 

Note that Pn(s) is optimal (and unique, of course) since the equation for 
constructing the PA imply that Pn{s) reproduces the first 2N moments of the 
distribution function, which we know. Using a diagonal PA of a Hamburger 
function to identify the optimal scales (ft) and weights (r*) is the basic idea 
behind the method of Gaussian quadrature for numerical integration . 

In the Brodsky-Lepage-Mackenzie (BLM) approach || one evaluates the 
coupling-constant at a scale that corresponds to the average momentum of 
the exchanged gluon (the BLM scale). This is exactly equivalent to approxi- 
mating the distribution function by a single 5-function located at its center. 
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This same effect can also be achieved simply by using a a; [0/1] PA for the 
leading (3q series [§J. In the method described above one uses an x[N — 1/N] 
PA of the leading (3q series that corresponds to approximating the momen- 
tum distribution function by a set of iV weighted 5 functions. Therefore it 
can be viewed as a generalization of the BLM method for higher-orders. 

In the following, I illustrate the above ideas with the vacuum-polarization 
D-function. The all-order distribution function p(s) [TIJ is plotted in fig. 6 as 
a continuous line. We see that there are contributions from both UV scales 
(positive s) and IR scales (negative s). The different symbols correspond 
to the locations (qi) and weights (r») of the diagonal PA's poles. For the 
x[0/l] PA, there is one ^-function at the BLM scale. For the a; [1/2] PA, 
there are two 5- functions, and so on. In Fig. 7 we see how using x[N — /N] 
PA corresponds to approximating the integral distribution function 0(s) by 
a piecewise constant function, composed of N steps. 

To conclude, we saw that diagonal PA's can be used to resum certain 
all-order effects that are related to the running of the coupling constant, 
and thus provide a systematic method for obtaining reliable scale-invariant 
predictions. There is a rigorous relation between diagonal PA and the mo- 
mentum distribution of virtual gluons. I stress that this result holds only 
for a single gluon exchange, i.e. within the large (3q approximation. The 
way to go beyond this approximation is still unclear. Nevertheless, from the 
results presented here for the Bjorken Sum-Rule, it is clear that PA's are an 
important tool for QCD phenomenology. 
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